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A turbulent-laminar banded pattern in plane Couette flow is studied numerically. This 
pattern is statistically steady, is oriented obliquely to the streamwise direction, and has 
a very large wavelength relative to the gap. The mean flow, averaged in time and in 
the homogeneous direction, is analysed. The flow in the quasi-laminar region is not 
the linear Couette profile, but results from a non-trivial balance between advection and 
diffusion. This force balance yields a first approximation to the relationship between 
the Reynolds number, angle, and wavelength of the pattern. Remarkably, the variation 
of the mean flow along the pattern wavevector is found to be almost exactly harmonic: 
the flow can be represented via only three cross-channel profiles as U(x, y, z) w XJo(y) + 
U c (y) cos(kz) + XJ s (y) sin(fcz). A model is formulated which relates the cross-channel 
profiles of the mean flow and of the Reynolds stress. Regimes computed for a full range 
of angle and Reynolds number in a tilted rectangular periodic computational domain 
are presented. Observations of regular turbulent-laminar patterns in other shear flows - 
Taylor-Couette, rotor-stator, and plane Poiseuille - are compared. 



1. Introduction 

Pattern formation is associated with the spontaneous breaking of spatial symmetry. 
Many of the most famous and well-studied examples of pattern formation come from fluid 
dynamics. Among these are the convection rolls which spontaneously form in a uniform 
layer of fluid heated from below and the Taylor cells which form between concentric 
rotating cylinders. In these cases continuous translational symmetries are broken by the 
cellular flows beyond critical values of the control parameter - the Rayleigh number or 
Taylor number. 

A fundamentally new type of pattern has bee n discovered in large-aspect-ratio shear 



flows in recent years by researchers at GIT-Saclay (Prigcnt fe Dauchotl . l2000tlPrigent et al 
I2002L l2003t IPrigent fe Dauchotl 120051 : iBottin et all Il998t ). Figure fTI shows an example 
from plane Couette experiments performed by these researchers. One sees a remarkable 
spatially-periodic pattern composed of distinct regions of turbulent and laminar flow. 
The pattern itself is essentially stationary. The pattern wavelength is large compared 
with the gap between the plates and its wavevector is oriented obliquely to the streamwise 
direction. 
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Figure 1. Photograph of a turbulent-laminar pattern in plane Couette flow from the Saclay 
experiment. Light regions correspond to turbulent flow and dark regions to laminar flow. The 
striped pattern of alternating laminar and turbulent flow forms with a wavevector k oblique 
to the streamwise direction. The wavelength is approximately 40 times the half-gap between 
the moving walls. The lateral dimensions are 770 by 340 half-gaps and the Reynolds number is 
Re = 385. Figure reproduced with permission from Prigent et al. 



The pattern emerges spontaneously from featureless turbulence as the Reynolds num- 
ber is decreased. This is illustrated in figure[5]with time series from our numerical simula- 
tions of plane Couette flow for decreasing Reynolds number (conventionally defined based 
on half the velocity difference between the plates and half the gap) . At Reynolds number 
500, the flow is uniformly turbulent. Following a decrease in the Reynolds number below 
400 (specifically 350 in figure [5]) the flow organises into three regions of relatively laminar 
flow and three regions of more strongly turbulent flow. While the fluid in the turbulent 
regions is very dynamic, the pattern is essentially steady. 

Shear flows exhibiting regular coexisting turbulent and laminar regions have a been 
know n for many years. In the mid 1960's, a state known as spiral turbulence was discov- 
ered (IColesL 119651 Ivan Attal . 1 19661 : IColes fc van Attal . Il966h in counter-rotating Taylor- 
Couette flow. Consisting of a turbulent and a laminar region, each with a spiral shape , 
spiral turbulence was further studie d in t he 1980s ( Andereck. Liu fc Swinnevl . Il986i 



Hegseth. A ndercck. Havot fc Pomeaul. 



- .19891) . Experiments by the Saclay re searchers 
(jPrigent fc Dauchotl . l2000l ; |Prigent et aLl . l200l 12003 : iPrigent fc Dauchotll2005l ) in a very 
large aspect-ratio Taylor-Couette system have shown that the turbulent and laminar re- 
gions in fact form a periodic pattern, of which the original observations of Coles and van 
Atta co mprised only one wavel ength. Analogues of these states occur in other shear flows 
as well. ICros fc Le Gall (|2002[ ) discovered large-scale turbulent spirals in the shear flow 
between a stationary and a rotating disk. Tsukahara. Seki. Kawamura fc Tochiol ( 20051 ) 
observed oblique turbulent-laminar bands in plane Poiseuille flow. A unified Reynolds 
number based on th e shear and the half-gap can be defined for these different flows 
( Prigent et al. , 120031) and is described in the Appendix. When converted to compara- 
ble quantities in this way, the Reynolds-number thresholds, wavelengths, and angles are 
similar for all of these turbulent patterned flows. The patterns are always found near the 
minimum Reynolds numbers for which turbulence can exist in the flow. 

In this paper we present a detailed analysis of these turbulent-laminar patterns. We 
will focus on a single case - the periodic pattern at Reynolds number 350. From computer 
simulations, we obtain the flow and identify the symmetries of the patterned state. We 
consider in detail the force balance responsible for maintaining the pattern. From the 
symmetries and harmonic content we are able to reduce the description to six ordinary- 
differential equations which very accurately describe the patterned mean flow. 
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Figure 2. Space-time plot from numerical simulations of plane Couette flow showing the 
spontaneous formation of a turbulent-laminar pattern at Re — 350. The kinetic energy in the 
mid-plane is sampled at 32 equally spaced points along an oblique cut (in the direction of pattern 
wavevector) through three wavelengths of the pattern. At time zero, Re = 500 and the flow is 
uniformly turbulent. Over about 3000 time units Re is decreased in steps to 350, and then held 
constant. 



2. Preliminaries 

2.1. Geometry 

The unusual but key feature of our study of turbulent-laminar patterns is the use of 
simulation domains aligned with the pattern wavevector and thus tilted relative to the 
streamwise-spanwise directions of the flow. Figure [3] illustrates this and defines our co- 
ordinate system. In figure [3fa) a simulation domain is shown as it would appear relative 
to an experiment, figure [TJ in which the streamwise direction (defined by the direction of 
plate motion) is horizontal. The near (upper) plate moves to the right and the far (lower 



plate t o the left in the figure. As we have discussed in detail (jBarklev fe Tuckerman . 



200530), simulating the flow in a tilted geometry has advantages in reducing computa- 



tional expense and in facilitating the study of pattern orientation and wavelength selec- 
tion. The important point for the present study is that the coordinates are aligned to 
the patterns. The z direction is parallel to the pattern wavevector while the x direction 
is perpendicular to the wavevector (compare figure [3Ja) with figure [1]). 

Figures [3jb) and (c) show the simulation domain as it will be oriented in this paper. In 
this orientation the streamwise direction is tilted at angle 9 (here 24°) to the x direction. 
This choice of angle is guide d by the experimental results a nd by our previous simula- 
tions. (In past publications ( Barklev fe Tuckerman . 2005 ai b ) we have used un-primed 



x — z coordinates for those aligned along spanwise-streamwise directions and primes for 
coordinates tilted with the simulation domain. Here we focus exclusively on coordinates 
fixed to the simulation domain and so for convenience denote them without primes.) In 
these tilted coordinates, the streamwise direction is 

e x cos 9 + e z sin 9 = ae x + /3e z (2.1) 

where 

a = cos6> = cos(24°) = 0.913, (3 = sin 9 = sin(24°) = 0.407. (2.2) 
We t ake L x = 10, fo r the reasons explained inljimenez fe Moinl ( 199l[ ); lHamilton et all 



(jl995l ); IWaleffd (|2003h : iBarklev fc Tuckerman! ((2005a|l6|) . Essentially, Z^sinfl must be 



near 4 in order to contain one pair of streaks or spanwise vortices, which are necessary to 
the maintenance of low Reynolds number wall-bounded turbulence. Although our simu- 
lations are in a three-dimensional domain, we will average the results in the homogeneous 
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Figure 3. Computational domain oriented at angle 6 to the streamwise-spanwise directions. 
The z direction is aligned to the pattern wavevector while the x direction is perpendicular to 
the pattern wavevector. The turbulent region is represented schematically by hatching, (a) 
Domain oriented with streamwise velocity horizontal, as in figure [T] (b) Domain oriented with 
z horizontal, as it will be represented in this paper. In (a), (b) the near (upper) plate moves in 
the streamwise direction; the far (lower) plate in the opposite direction, (c) View between the 
plates. 



x direction, as will be explained in section l2~3l For most purposes it is sufficient to view 
the flow in the z — y coordinates illustrated in figure [3Jc). The midplane between the 
plates corresponds to y = 0. 

The length L z of our computational domain is guided by the experimental results and 
by our previous simulations. One of the distinctive features of the turbulent- laminar 
patterns is their long wavelength relative to the gap between the plates. A standard 
choice for length units in plane Couette flow is the half-gap between the plates. In the 
simulation with L z — 120 and 9 = 24° shown in figure [U a pattern of wavelength 40 
emerged spontaneously from uniform turbulence when the Reynolds number was lowered 
to Re — 350. For this reason, the simulations we will describe below are conducted with 
L z — = 40. The corresponding wavenumber is 

fc = ^ =0.157. (2.3) 
40 

This large wavelength, or small wavenumber, expresses the fact that the pattern wave- 
length in z is far greater than the cross-channel dimension. 

2.2. Equations and Numerics 
The flow is governed by the incompressible Navier-Stokes equations 

c^u 1 

= -(u • V)u - Vp+ — V 2 u infi, (2.4a) 
at Re 

= V • u in Q., (2.46) 

where u(x, t) is the velocity field and p(x, t) is the static pressure. Without loss of 
generality the density is taken to be one. The equations have been nondimensionalized 
by the plate speed and the half gap between the plates. O is the tilted computational 
domain discussed in the previous section. 

No-slip boundary conditions are imposed at the plates and periodic boundary condi- 
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tions are imposed in the lateral directions. In our coordinates the conditions are 



u(x, y — ±1, z) = ±(e x cos 9 + e z sin( 
u(x + L x ,y, z) = u(a;,y, z) 
u(x,y,z + Lg) = u(x,y,z). 



(2.5a) 
(2.56) 
(2.5c) 



Linear Couette flow u L is a solution to (|2.4|) — (|2.5p . which is stable for all Re and 
satisfies 



V 2 u L = (u L • V)u L = 

In our tilted coordinate system, 

u L = y(e x cos 9 + e z sin0) — y(ae x + /3e z ) = u L e x + w L e 2 



(2.6) 



(2.7) 



The Navier-Stokes equations (|2.4[) with boundary c ondition s (12.51) are simulated usin g 
the spectral-element (x-y) - Fourier (z) code Prism ( Henderson fe Karniadakisl 1995 ). 
We use a spa tial resolution consistent with previous studies ((Hamilton et al 1 119951 : 



Waleffd . 120031 ). Specifically, for a domain with dimensions L x = 10 and L y — 2, we 



use a computational grid with 10 elements in the x direction and 5 elements in the y 
direction. Within each element, we use 8th order polynomial expansions for the primitive 
variables. In the z direction, a Fourier representation is used and the code is parallelized 
over the Fourier modes. Our domain with L z = 40 is discretized with 512 Fourier modes 
or gridpoints. Thus the total spatial resolution we use for the L x x L y x L z = 10 x 2 x 40 
domain can be expressed as N x x N y xN z = 81 x 41 x 512 = 1.7 x 10 6 modes or gridpoints. 



2.3. Dataset and averaging 

The focus of this paper is the mean field calculated from the simulation illustrated by the 
spatio-temporal diagram in figure |4ja). The velocity field in the portion of the domain 
shows high-frequency and high- amplitude fluctuations, while the flow in the right portion 
is basically quiescent. We will call the flow on the left turbulent, even though it could 
be argued that it is not fully developped turbulence. We will call the flow on the right 
laminar, even though occasional small fluctuations can be seen in this region. 

The turbulent-laminar pattern subsists during the entire simulation of 14 x 10 3 time 
units. However the pattern undergoes short-scale "jiggling", seen particularly at the 
edges of the turbulent regions, and longer-scale drifting or wandering in the periodic 
z direction. We seek to describe the field which results from smoothing the turbulent 
fluctuations, but for which drifting is minimal, by averaging over an appropriate time 
interval. The desired averaging time interval represents a compromise between the short 
and long timescales. We have chosen to average the flow in figure @|a) over the shaded 
time interval [t, t + T] = [6000, 8000] , during which the pattern is approximately station- 
ary. 

The time-averaged flow is homogeneous in ^-direction. This is illustrated in figure[4jb) 
where we plot one of the velocity components time-averaged flow over the interval 
[6000,8000]. Cuts at different x locations show that there is essentially no variation 
in the a;-direction. All other quantities are similarly independent of x. It is therefore 
appropriate to consider mean flows as averages over the x direction as well as over the 
time. 
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We define mean flows as 

rt+T r Lx 



1 1 r +J [ u » 

(u)(y,z) = ——J J u(x,y,z,t) dxdt (2.8a) 



^ 1 rt+T pL x 



{p)(y yZ ) = -— / p(x,y,z,t)dxdt. (2.86) 

1 i^x Jt JO 

The mean fields obey the averaged Navier-Stokes equations 

= - «u) • V) (u) - <(u • V) u) - V(p> + -^V 2 (u) (2.9a) 

tie 

= V-(u), (2.96) 

where 

u = u (u) (2.10) 
is the fluctuating field and () denotes x-t average. The mean fields are subject to the 



same boundary conditions as equations (|2.4j) . We denote the Reynolds-stress force from 
the fluctuating field in equations (|2.9I) by F: 



F = -((u • V) u) = -V • (uu), (2.11) 

We shall focus almost exclusively on the difference between the mean flow and linear 
Couette flow, for which we introduce the notation 

U=(u)-u L , (2.12) 

as well as P = (p). 

Letting the components of U be denoted by (U, V, W) and the components of F be 
denoted by (F U ,F V ,F W ), then the averaged Navier-Stokes equations for the deviation 
from linear Couette flow in component form become 

= - (Vd y + (W + /3y) d z )(U + ay) + j^{d 2 y + 8 2 Z )U + F u (2.13a) 

= - (Vd y + (W + f3y) d z )V - d y P + -^(d 2 y + d 2 z )V + F v (2.136) 

= - (Vd y + (W + 0y) d z ) [W + [3y) - d z P + -±-(6$ + d 2 z )W + F w (2.13c) 

= d y V + d x W. (2.13d) 

U is required to satisfy homogeneous boundary conditions at the plates 

U(y = ±l,z) = (2.14) 

and periodic boundary conditions in z. 

A system of this type, with three components depending on two coordinates, is some- 
times called 2.5 dimensional. The transverse, or out-of-plane flow U(y, z) appears only in 
the first equation and is effectively a passive scalar advected by the in-plane flow (V, W) 
and driven by the Reynolds-stress force F u . The in-plane flow can be expressed in terms 
of a streamfunction ^ where 

Ve y + We z = e x x V* = -d z fe y + d y fe z . (2.15) 



We shall use both (U, V, W)(y, z) and (U, ^ / )(y, z) to describe the mean flows. 
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Figure 4. (a) Timeseries of a turbulent-laminar pattern. Shown is the kinetic energy E = u-u/2 
along the line x — y — at 32 equally spaced points in z for ^ t ^ 14000. The interval 
[6000,8000] used for time averaging is shown in grey, (b) Time-averaged velocity at five x 
locations illustrating the ^-independence of the time-averaged flow. Plotted is ^ J dt(u — u L ), 
the average x component of velocity with linear Couette flow subtracted, averaged over the 
interval [6000,8000] indicated in (a). Color range from blue to red: [-0.4, 0.4]. 




Figure 5. U(y,z): transverse component of mean flow, ^{y, z): streamfunction of in-plane 
mean flow. A long cell extends from one laminar-turbulent boundary to the other. Gradients 
of ^ are much larger in y than in z, i.e. \W\ ^> |V|. In the laminar region at the center, 
W, V ~ 0. E t uib(y, z): mean turbulent kinetic energy (u • u)/2. There is a phase difference of 
A z /4 = 10 between extrema of -Eturb and of U. P(y, z): mean pressure field. Pressure gradients 
are primarily in the y direction and within the turbulent region. Color ranges for each field from 
blue to red: U [-0.4, 0.4], * [0, 0.09], E tulh [0, 0.4], P [0, 0.007] . 

3. Results 

We present a characterisation of the turbulent-laminar pattern at Re = 350. We 
describe in detail the mean flow, its symmetries, and the dominant force balances within 
the flow. Our goal here is not to consider closures for averaged Navier-Stokes equations 
(|2.13| . We will make no attempt to model the turbulence, i.e. to relate the Reynolds- 
stress tensor (uu) to the mean flow U. Instead we use fully resolved (three-dimensional, 
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timc-dcpcndcnt) numerical simulations of the turbulent flow to measure both the mean 
field U and Reynolds-stress force F. From these we extract the structure of these fields 
and the dominant force balances at play in sustaining turbulent-laminar patterns. 

3.1. Mean flow 

The mean flow is visualised in figure [5] via the transverse, out-of-plane flow U(y,z) and 
the in-plane streamfunction ^(y,z). Recall [equation (|2.12p ] that these fields are the 
deviations of the mean flow from linear Couette flow u L . The mean turbulent kinetic 
energy 

£ turb = -(u-u) (3.1) 

serves to clearly identify the turbulent region. In these and subsequent plots, the middle 
of the laminar region is positioned at the centre of the figure and the turbulent region at 
the periodic boundaries of the computational domain. In figure [5] (but not in subsequent 
figures), plots are extended in the z-direction one quarter-period beyond each periodic 
boundary to help visualise the flow in the turbulent region. The pattern wavelength is 
X z = 40, so that z = 30 and z = — 10 describe the same point, as do z = —30 and z = 10. 

The mean flow can be described as follows. U is strongest in the turbulent-laminar 
transition regions. In the transition region to the left of centre (z = —10) in figure [5) 
U is negative and primarily in the upper half of the channel. To the right of centre 
(z = 10), U is positive and is seen primarily in the lower half of the channel. Comparison 
with turbulent kinetic energy shows that the transverse mean flow U is out of phase 
with respect to the fluc tuating field u by This has be en seen experimentally by 

Coles fc van Attal (|l966h and Prigent et al. (|200ll2003ll2005h . 



The in-plane flow \P in figure [S] has a large-aspect ratio cellular structure consisting of 
alternating elliptical and hyperbolic points. The flow around the elliptical points, located 
in the centre of the turbulent regions, rotates in a counter-clockwise sense, opposing linear 
Couette flow. In the vicinity of the hyperbolic points, centred in the laminar regions, the 
in-plane deviation from linear Couette flow is very weak (W and V nearly zero). 

Figure [S] shows y-profiles at four key points equally spaced along the pattern: centre 
of the laminar region, turbulent-laminar transition region, centre of the turbulent region, 
and the other turbulent-laminar transition region. While the V profile is plotted, its 
variation is very small on the scale of U and W and can essentially be used to indicate 
the axis. Figure [7] shows profiles for the full mean flow (u) = U + u L containing the 
linear Couette profile. 

The U profiles in figure [7] are S-shaped, of the type found in turbulent Couette flow. 
This is to be expected in the turbulent region, even at these low Reynolds numbers. 
However, it is very surprising that the U profile in the laminar region is also of this form. 
In the laminar region, local Reynolds stresses are absent (see figure [5]) and so cannot 
be responsible for maintaining the S-shaped velocity profile in the laminar regions. The 
other prominent features in figures [6] and [7] are the asymmetric profiles at the transition 
regions. 

The relationship between the mean flow field and the regions of turbulence can be seen 
in figure [8] Here the flow is shown in the standard orientations. In each view, greyscale 
indicates the size of the turbulent energy and the arrows show the mean flow within the 
plane. In the top two views, the flow is shown in the streamwise-spanwise planes located 
at y = 0.725 and at y = —0.725. The next view shows the flow between the plates, i.e. in a 
streamwise-cross-channel plane, and the last shows an enlargement of one of the laminar- 
turbulent transition regions. Note that the length L z = 40 of our tilted computational 
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Figure 6. Mean flow profiles in y at four equally spaced locations in z. From left to right: centre 
of the laminar region (z = 0), laminar-turbulent boundary (z = 10), centre of the turbulent 
region (z — 20) and turbulent-laminar boundary (z = —10). Components U (blue, solid), V 
(green, dotted), W (red, dashed) of deviation from linear Couette flow u L . In the laminar 
region, W ~ 0, indicating no deviation from u L . V is very small throughout. 
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Figure 7. Same as in figure [6] but with laminar Couette flow u L included. 



domain corresponds to a streamwise length of L z j sin 9 = 40/. 407 = 98.3 « 100 and to a 
spanwise length of L z / cosO = 43.78. 

The flow in figure[5]can be compared with the mean flow reported bv lColes &: van Atta 
(1966) in experiments on turbulent spirals in Taylor-Coutte flow. Coles and van Atta 
measured the mean flow near the midgap between the rotating cylinders and noted an 
asymmetry between the mean flow into and out of turbulent regions. They found that 
the mean flow into turbulent regions was almost perpendicular to the turbulent-laminar 
interface whereas flow out of the turbulent region was almost parallel to the turbulent- 
laminar interface. We also observe a striking asymmetry between the mean flow into and 
out of the turbulent regions. The orientation of our mean flow does not agree in detail 
with that of Coles and van Atta, but this is most likely due to the fact that Coles and van 
Atta considered circular Taylor-Couette flow and measured the flow near the mid-gap. 
Referring to figures [5] and [6] one sees that the mid-plane (y = 0) is not the ideal plane 
on which to obverve the mean flow since its structure is most pronounced between the 
midplane and the upper or lower walls. 

Before considering the symmetries and force balances in detail, it is instructive to 
consider the dominant force balance just at the centre of the laminar region. Recall that 
one of the more interesting features of the mean flow is that the U profile appears very 
similar to a turbulent profile, even in the absence of turbulence in the laminar region. 
Here the balance is dominated by advection and viscous diffusion, as shown in figure 03 
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Figure 8. Mean velocity components seen in three planes with standard orientation for Couette 
flow. The turbulent regions are shaded. Top: velocity components in the streamwise-spanwisc 
plane at y = 0.725 (upper part of the channel). Middle: same except y — —0.725 (lower part 
of the channel). Bottom: flow in a constant spanwise cut. The mean velocity is shown in the 
enlarged region. 



Equation (|2.13n[) for flow in the x-direction is 

= - (Vd y + (W + (3y)d z ) (U + ay) + ^-(^ + d 2 z )U + F u . (3.2) 

Variations in y dominate variat ions in z, i.e . the usual boundary-layer approximation 
{&l + d 2 z )U ~ d 2 U holds; see, e.g. IPoDel (l2000h . Indeed, approximating the y dependence 
of U by the functional form sin(7ry) suggested by figure El we have 



O 



(2V40) 2 



400 



(3.3) 



This is confirmed by the second panel of figure [9] In the centre of the laminar region F u , 
V, and W are all negligible, so that —(3yd z U dominates the advective terms, as shown 
in the third panel of figure [H Thus the balance between advection and viscosity in the 
laminar region is 

Pyd z U^^d 2 y U. (3.4) 
This equation is appealingly simple and yet leads immediately to some interesting 
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Figure 9. Balance of forces in the center of the laminar region. Left: Forces in the U direc- 
tion. Advective terms (blue, solid), viscous terms (red, dashed), Reynolds-stress terms (green, 
dotted). Middle: Viscous terms in the U direction. (1/ Re) dyU (red, dashed) dominates dom- 
inates (l/Re)d z U ( green, dotted). Right: Advective terms in the U direction. Curves show 
—f5yd z U (blue, solid) —Wd z U (black, dash-dot), — aV (green, dotted), —Vd y U (red, dashed). 
Right-most: Advective terms in the W direction (for later reference). Curves show — f3yd z W 
(blue, solid) -Wd z W (black, dash-dot), -/3V (green, dotted), -Vd y W (red, dashed). 



conclusions. The first is that a non-zero tilt angle 9 is necessary to maintain the S- 
shaped U profile in the laminar region, since otherwise [3 = sin# = and U could be at 
most linear in y and would in fact be zero, due to the homogeneous boundary conditions 
(|2.14p . The second conclusion follows from consideration of y parity. The multiplication 
by y on the left-hand-side reverses y-parity, while the second derivative operator on the 
right-hand-side preserves y parity. The conclusion is that U should be decomposed into 
odd and even components in y and equation (|3.4[) should actually be understood as two 
equations coupling the two components. Specifically, as can be seen in figure [HI U is odd 
in y in the centre of the laminar region, yet d z U must be even for equation (|3.4p to hold. 

The remainder of the paper is devoted to formalising, demonstrating and extending 
this basic idea. 

3.2. Symmetry and Fourier modes 

We now consider in depth the symmetry properties of the flow. We start with the 
symmetries of the system before averaging, that is, the Navier-Stokes equations (|2.4|) 
and boundary conditions (|2.5|) . The system has translation symmetry in x and z as well 
as centrosymmetry under combined reflection in x, y and z: 

K xyz (u, v, w)(x, y, z) = (-u, -v, -w)(x - x, -y, z Q - z) (3.5) 

where the origin xo, zq can be chosen arbitrarily. Linear Couette flow u L possesses all 
the system symmetries, as does the mean flow at Reynolds numbers for which turbulence 
is statistically homogeneous in x and z. 

Note that in the absence of tilt (9 = 0), the system possesses two reflection symmetries. 
These can be taken to be n xyz and reflection in the spanwise direction. For the tilted 
domain (at angles other than multiples of 90°), the only reflection symmetry is n xyz . This 
can be seen in figure [3Ja): for general tilt angles 9, spanwise reflection does not preserve 
the domain, i.e. does not leave the periodic boundaries in place. The experimental 
system shown in figure Q] possesses spanwise reflection symmetry and hence bands can be 
observed in the either of the two symmetrically related angles, the choice is dictated by 
factor such as initial conditions. By design, our tiled computational domain precludes 
the symmetry-related pattern given by spanwise reflection. 

The transition to the turbulent-laminar patterned state breaks symmetry. Specifically, 
both the mean flow and the Reynolds-stress force break z-translation symmetry but break 
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z even 


z odd 


y even 


0.03% 


25.48% 


y odd 


74.48% 


0.01% 



Table 1. Energy J dx J dy J dz |U| 2 /2 of deviation from Couette flow contained in modes 
with different symmetries in y and z. Modes with centrosymmetry (opposite parity in y and z) 
contain 74.48%+25.48%=99.96% of the total energy. Reflection in z is about the centre of the 
laminar region. 



neither x-translation symmetry nor centrosymmetry. The spatial phase of the pattern in 
z is arbitrary, but given a phase there are two values of zq, separated by half a period, 
for which the flow is invariant under n xyz , as is typical for a circle pitchfork bifurcation 
( Crawford fc Knobloch . 199lh . As can be seen in figure 03 the values of zq about which 
the patterns are centrosymmetric are the centres of the laminar (zq = 0) and of the 
turbulent (z — ±20) regions. 

The centrosymmetry operator for our averaged fields U, which depend only on y and 
z, is 

Kyz (U, V, W)(y, z) = {-U, -V, -W)(-y, z - z) (3.6) 

Since the Reynolds-stress force (F u ,F V , F w ) is centrosymmetric in the case we consider, 
then the averaged equations (|2. 13|) for the mean field have centrosymmetry. 

We formalise this further as follows. Any ^-independent field g can be decomposed 
into even and odd functions of y and z as 

g(y, z) = g ao (y, z) + g oe (y, z) + g eo (y, z) + g ee (y, z) (3.7) 

where, for example, g oe is odd in y and even in z — Zq. Applying the operator in (|3 . 6|) to 
(|3.7p . we obtain 

K yzg{y, z) = -goo(-y, Zq-z)- goe(-y, Z Q -Z)- geo(-y, Zq-z)- gee{~y, Z Q - z) 

= -gooiV, z) + goe{y, Z) + geo(y, z) - gee(y, z) (3.8) 

For the field g to be centrosymmetric requires n yz g — g, so that in fact 

g{y,z) = g oe (y,z) + g eo (y,z) (3.9) 

Table [I] as well as figure shows that this is indeed the case for U; it holds for F as 
well. 

We now Fourier transform in z to further decompose the mean velocity and the 
Reynolds-stress force. We find that the z-wavenumbers and ±fc have contributions 
to U which are an order of magnitude higher than the remaining wavenumber combina- 
tions. See tabled The deviation from the z average is thus almost exactly trigonometric, 
with almost no higher harmonic content. The dominance of these terms in the Fourier 
series means that U and F can be represented by only three functions of y, namely: 

g(x, y, z) = g (y) + g c (y) cos(fcz) + g s (y) sin(fcz) (3.10) 

which is a special case of (|3.9[> . with the first two terms of (|3.10p coinciding with g oe {y, z) 
and the last to g eo . Thus, go and g c are odd functions of y, while g s is even. The 
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z wavenumber 


z even (cosine) 
k > 2fc 


2 odd (sine) 
k ^2k 


y even 






25.2% 


0.3% 


y odd 




69.7% 




4.7% 


0.1% 





Table 2. Energy contained in 2 Fourier modes. Modes retained (in boxes) are Uo(y), 
U c (y)cos(kz) and U s (j/) sin(fcz). These contain 69.7%+4.7%+25.2%=99.6% of the total 

energy. 



cos(0z)=l cos(kz) sin(kz) 




-0.2 0.2 -0.2 0.2 



Figure 10. Fourier decomposition of mean velocity. U component (blue, solid), V component 
(green, dotted), W component (red, dashed). W c « —Wo, corresponding to the fact that W 
shows no deviation from the linear in the laminar region. 



fields thus consist of a z-independent component go and two components which vary 
trigonometrically and out of phase with one another, g c dominating in the laminar and 
turbulent regions and g s dominating in the boundaries between them. Moreover, g s 
dominates in the bulk, since go and g c are odd in y and thus zero in the channel centre. 



g = So(y) + Sc{y) z = : Centre of laminar region (3.11a) 

g = So(y) + Ss(y) z = A z /4 = 10 : Laminar-turbulent boundary (3.116) 

g = go(y) — gc(y) z = Ax/2 = 20 : Centre of turbulent region (3.11c) 

g = So(y) — Ss(y) z — 3A 2 /4 = 30 : Turbulent-laminar boundary (3. lid) 

Figure [TU] shows the three trigonometric components, each a function of y, obtained 
by Fourier transforming U,V, and W. Figure QT] shows U, V and W as functions of z at 
locations in the upper and lower channel and compares them with the values obtained 
from the trigonometric formula (|3 . 10[) using the functions shown in figure [TU1 Figures [T2l 
[TBI and [TU depict U(y,z), ty(y,z) and F u (y,z) with their trigonometric decompositions. 
Each of these figures uses only the three scalar functions of y, figure [TD1 to reproduce the 
corresponding two-dimensional field. As shown by equation (|2.15p . the streamfunction 
4" of a centrosymmetric field has symmetry opposite to that of the velocity components, 
i.e. it is composed of functions of the same parity in y and z. 

Figures [J5] and QJj] show the three Reynolds-stress forces and their Fourier decompo- 
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Figure 11. Mean flow as a function of z at y = 0.725 (lower curves) and y = —0.725 (upper 
curves). U (blue, solid), V (green, dotted), W (red, dashed). Dots show values calculated from 
trigonometric formula (|3.10[) . 



sitions. Each component obeys F c a — Fo, a necessary condition for F to vanish at 
the centre of the laminar region, as shown by equation (|3.11al) and also illustrated in 
figure [HI More precisely, 



F l c J = -1.09^, 



F 



-1.22F y , = -1.16^ 



In addition, 

F a -9 y (uv). 

as is typical for turbulent channel flows; see, e.g. Pope (2000). 



(3.12) 
(3.13) 



3.3. Force balance for U 

We now turn to understanding the balance of forces responsible for maintaining the 
mean flow profiles. We focus primarily on U, both because it is the component of largest 
amplitude and also because it appears only in equation (|2.13aj) : U is subject to Reynolds- 
stress and viscous forces, and is advected by (V, W) but is not self-advected. We begin 
by showing the balance of forces in the U direction as a function of z at locations in the 
upper and lower channel in figure 1171 One can again see the centrosymmetry of each 
of the forces, i.e. invariance under the combined operations of reflection in y and z and 
change of sign. The Reynolds-stress force disappears at the center of the laminar region 
and the advective and viscous forces exactly counterbalance, as emphasized in the figures 
on the right. Figure IT51 shows another view of this balance, displaying the forces as a 
function of y at four locations in z. As previously stated, V 2 C/ is dominated by dyU and 
F u by —d y (uv). In figure [T9l we show the Fourier-space analogue of figure [T8l 

We now turn to the more complex advective forces, whose Fourier decompositions are 
shown in figure [20l The cos(Oz) component of the advective force is small but non-zero. 
Because this term results from the product of trigonometric functions, it also provides 
a measure of the generation of higher harmonics, a point which we will explore further 
in section 13.51 The advective cos(fcz) term is well approximated by the contribution 
from advection by w L — j3y. The advective sin(fcz) term is dominated near the walls by 
advection by w L , but in the bulk by advection by V. Properties of the cos(fcz) and sin(fcz) 
modes echo their physical space counterparts: the advective term is well approximated 
by advection by w L in the laminar region, as was shown in figure O while the advective 
forces in the laminar-turbulent boundaries combine advection by w L near the walls and 
by V in the bulk. 

We illustrate these conclusions via schematic visualizations of the dynamics of U. Fig- 
ure[2T]illustrates the dynamics in the laminar and turbulent regions. The dynamics in the 




Figure 12. Mean velocity U and its trigonometric decomposition. Because the magnitude of 
the fields vary, different color scales are chosen to emphasize qualitative features. U (—0.2, 0.2), 
U (-0.2,0.2), Uc (-0.05,0.05), U a (-0.154,0.154). 
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Figure 13. Mean streamfunction ^(y, z) for deviation of in-plane flow from linear Couette flow 
and its trigonometric decomposition. Color scale is * (0, 0.09), *o (0, 0.046), * c (-0.042, 0.042), 
* s (-0.008,0.008). 




+ 
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2:; fo o iq so 

Figure 14. Reynolds-stress force F u and its trigonometric decomposition. Color scale is 
F u (-0.017,0.017), Fo (-0.0085, 0.0085), F? (-0.0085, 0.0085), ^(-0.0085,0.0085). 



laminar region are essentially described by the simple balance between viscous diffusion 
of U profiles and advection by linear Couette flow in z, given by equation (|3.4[) . Viscous 
diffusion tends to reduce curvature, but the profiles have greater curvature upstream 
(to the left for the upper channel, to the right for the lower channel). Hence advection 
replenishes the curvature damped by viscosity. However, this trend towards greater cur- 
vature upstream cannot continue indefinitely, since the pattern is periodic in z. Hence 
eventually a maximum is reached (at a turbulent-laminar boundary), beyond which the 
curvature decreases upstream. Thus, in the turbulent region, advection and diffusion act 
together to decrease curvature and must both be counter-balanced by turbulent forcing. 
These features are essentially described by the cos(Oz) = 1 and cos(fcz) modes. Figure l22l 
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Figure 15. Reynolds-stress force F = — ((u • V)u) as a function of y. Curves show F u (blue, 
solid), F v ( green, dotted) and F w (red, dashed). 
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Figure 16. Fourier decomposition of Reynolds-stress force. F u components (blue, solid), F v 
components (green, dotted), F w components (red, dashed). F c ~ — Fo, as required for the 
vanishing of F at z — 0. 



illustrates the dynamics in the turbulent-laminar boundaries. These dynamics include 
advection by V in the bulk, leading to the U > (U < 0) patch in the lower right (upper 
left) of figure fT2l and are described by the sin(fcz) mode. 

3.4. Force balance for W and V 

Figure [23] shows the balance of forces in the W direction and figure EH its analogue in 
Fourier space. This balance resembles that in the U direction shown in figures fT8l and [T9l 
In physical space (compare the leftmost panels of figures |2"31 and [TS]) . the main difference 
is that the advective and viscous forces are both small in the laminar region, in keeping 
with the fact that W ~ 0. The pressure gradient d z P is far smaller than the other forces 
throughout (see below). In Fourier space (compare the middle panels of figures l24l and 
fl9|) , the main difference with the U balance is that the relative importance of the viscous 
and advective forces in the cos(fcz) balance is reversed from that in the case of U : for 
W, the viscous component is larger than the advective component, which is especially 
small in the bulk. The decomposition of the advective terms (figure [25]) shows that, as is 
the case for U, the advective cos(fcz) term is well approximated by the contribution from 
advection by w L = f3y, whereas all four advective components contribute to the sin(fcz) 
term. 

The balance of forces in the V direction is entirely different. The dominant balance in 
this equation is: 

= -d y P + F v (3.14) 



as shown in figure 1261 This is typical for turbulent channel flows; see, e.g. (jPopel . 1200(1 




Figure 17. Mean forces in U direction as a function of z at y = ±0.725 for turbulent-laminar 
pattern at Re = 350. Advective — (U • V)t/ (blue, solid), viscous V 2 C7 (red, dashed), and 
turbulent — ((u ■ V)w) (green, dotted) forces. In the laminar region (z m 0), the Reynolds-stress 
force vanishes and the viscous and advective forces are equal and opposite to one another. In 
figures on right, enlarged around the laminar region, \7 2 U and +(U- V)U are shown to emphasize 
equality between viscous and advective forces. 

This balance between the mean pressure gradient P and the Reynolds-stress force F v 
does not constrain or provide information about any of the velocity components. Since 

F v = -V • (ire) « ~d y (v 2 ) (3.15) 

we in fact have 

P«-(w 2 ) (3.16) 

up to a small z-dependent correction. Figure [5] shows the pressure field P calculated from 
(|3.16p and suggests that its y dependence can be approximated by the functional form 
cos(7ry/2). This leads to an estimate of the relative importance of the pressure gradients 
in the y and z directions: 



d y P\ _ it/2 



while our data shows 



°<rfp =1r* 10 (3 - 17) 



(d y P) max _ 0.012 



W. - 0.0017 - 7 -» 5 I 3 ' 18 ' 
The same estimate applies to the relative magnitudes of V and W, using the stream- 
function shown in figure [5] 



0\ — \ =0\ = —!— = 1Q (3.19) 

1 V ) \d z y ) 2tt/40 v ' 



while the actual ratio of maximum values is 

Wmax 0.15 



11. (3.20) 

Knax 0.013 V ^ 
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Figure 18. Balance of forces in the U direction. Curves show advective force (blue, solid), 
viscous force (red, dashed) and Reynolds-stress force (green, dotted). In the laminar region, 
the Reynolds-stress force is negligible and the advective and viscous forces counter-balance one 
another. 
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Figure 19. Balance of forces in the U direction, decomposed into modes. Curves show ad- 
vective (blue, solid), viscous (red, dashed) and Reynolds-stress (green, dotted) forces. Mode 1: 
Reynolds-stress and viscous forces approximately counterbalance each other, cos(fcz): advection 
is larger than viscous force, which is especially small in the bulk. 
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Figure 20. Advective terms in the U direction, decomposed into modes. Curves show —(3yd z U 
(black, solid) —Wd z U (blue, dash-dot), —aV (green, dotted), —Vd y U (red, dashed). The 1 
mode is generated by the product — kW c U s /2; a second harmonic of the same small size is also 
generated. The cos(fcz) mode is dominated by —/3ykU s - The sin(kz) term is dominated by 
fiykUc near the boundaries and by —V s d y (Uo + ay) in the bulk. 



3.5. Model equations 

We now derive a system of ordinary differential equations by substituting the trigonomet- 
ric form (|3.10p into the Reynolds-averaged Navier-Stokes equations (|2.13|) . The drawback 
in this procedure is the usual one, namely that this form is not preserved by multiplica- 
tion. However, Table [2] shows that higher harmonics contribute very little to U. 
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Figure 21. Schematic depiction of the dynamics of U near the centres of the laminar and of 
the turbulent regions. The cross-channel direction is exaggerated. Shown are U (profiles) and 
W + f3y (arrows). Viscous diffusion tends to diminish both peaks of the profile. In laminar region 
surrounding z — 0, the peaks in the upper half-channel increase in amplitude with decreasing z; 
advection towards positive z (upper arrow) replenishes these peaks, maintaining U. Conversely, 
the peaks in the lower half-channel increase with z; advection towards negative z (lower arrow) 
replenishes these peaks. That is, the sign of — (W + fly) d z U is opposite to that of dyU in both 
the upper and lower parts of the laminar region. In the turbulent region around z — ±20, 
the size of the upper (lower) peak decreases to the left (right) and so advection, like viscous 
diffusion, acts instead to diminish the peaks. U is maintained by the Reynolds-stress force, which 
counterbalances both. The effect is to modulate the amplitude of the U profiles periodically in 
z. 
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Figure 22. Schematic depiction of the dynamics of U near the turbulent-laminar boundaries. 
The cross-channel direction is exaggerated. Shown are U + ay (profiles), and (V, W + fly) 
(arrows). Near the upper and lower walls, the U + ay profiles are advected towards increas- 
ing/decreasing z by W + (3y. In the bulk, advection by V is significant. At z ~ 10, V advects 
downwards the right-moving fluid in the upper portion of the channel. At z ~ —10, V advects 
upwards the left-moving fluid in the lower portion of the channel. The effect is to tilt the {7 = 
boundary periodically in z. 

We expand the advective term as: 

((U + u L ) -V) (U + u L ) = (V d y + (W + 0y) d z )(V + aye x + /3ye z ) 

= (V c cos(fcz) + V s sin(fcz))(Uo + U' c cos(fcz) + U' s sin(fcz) + ae x + j3e z ) 
+ (W + (3y + W C cos(fcz) + W s sin(fcz))(-fc U c s'm(kz) + k U s cos(fc2:)) 



= - (V c U' c + V s U' s + k (W c U s - W s U c )) (3.21a) 

+ (V C (U' + ae x + /3e z ) + k (W + f3y)V s ) cos(fcz) (3.216) 

+ (V S (U' + ae x + (3e z ) - k (W + /3y)U c ) sin(fcz) (3.21c) 

+ 1 (V c V' c - V s U' s + k (W c U s + W s U c )) cos(2fcz), (3.21d) 



where primes denote y differentiation. We neglect the second harmonic term (|3.21d|) . 
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Figure 23. Balance of forces in the W direction. Curves show advective term (blue, solid), 
viscous force (red, dashed) and Reynolds-stress force (green, dotted). In the laminar region, 
W ~ and each of the forces is negligible. In the turbulent region, the viscous and Reynold- 
s-stress forces counter-balance one another. In the laminar-turbulent boundaries, the advective, 
viscous and Reynolds-stress forces all play a role. 
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Figure 24. Balance offerees in the W direction, decomposed into modes. Curves show 
advective term (blue, solid), viscous term (red, dashed) and turbulent forcing term (green, 

dotted). 
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Figure 25. Advective terms in the W direction, decomposed into modes. Curves show —f3yd z W 
(black, solid) -Wd z W (blue, dash-dot), -/3V (green, dotted), -Vd y W (red, dashed). The 
cos(fcz) mode is dominated by —f3ykW s . The sin(fcz) term is dominated by /3ykW c near the 
boundaries and by —V s d y (Wo + /3y) in the bulk. 



and will discuss the accuracy of this approximation below. We now rewrite the U and 
W components of the averaged momentum equations, neglecting the z-derivatives d^U, 
d\W and d z P, as justified by equations (|3.3p and (|3.17[) : 

= - (Vd y + (W + 0y) d z ) (U + ay) + ±-6%U + F u (3.22a) 

= - (Vd y + (W + 0y) d z ) [W + 0y) + ^d 2 y W + F w (3.226) 
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Figure 26. Forces in the V direction. Reynolds-stress force F (red, dashed) is counterbalanced 
by pressure gradient —d y P (blue, solid). Both are zero in the laminar region. Advective and 
viscous forces (green, dotted) are negligible throughout. 



Substituting (|3.21aj) — (|3.21c[) in (|3.22|) and separating terms in cos(Oz) = 1, cos(fcz) and 
sin(fcz), we obtain 



= ™ [V C U' C + V S U' S + k(W c U s - W S U C )] + -^f/o + Fq 



= -V c {% + a) - k (Wo + 0y)U s + -^-U'J 

tie 

= -V s (U' a + a)+k(W + f3y)U c + -±-U'J 

He 



F„ 



Ft 



w 



= - 5 [VcW^ + V S W' S ] + — < + F K 

= -V C (W^ + 0) - k (Wo + 0y)W s + -i-W? + FY 

He 

= -V S (W^ +0) + k(W o + 0y)W c + -1-W'J + Ff 

He 



(3.23a) 
(3.236) 
(3.23c) 
(3.23d) 
(3.23e) 
(3.23/) 



where the Fourier modes of V and W are related via those of the streamfunction ^ of 

V =0 Wo = % (3.24a) 

V e = VF C - ^ (3.246) 

V s = fc* c W s = *i (3.24c) 



and where homogeneous boundary conditions are imposed 

= ?7 = £7 C = U s &ty = ±1 
= W = W c = VF S at y 



±1 



(3.25a) 
(3.256) 



System (|3.23[) with boundary conditions (|3.25p is composed of six ordinary differential 
equations coupling the six scalar functions Uq, U c , U s , 'J'oj of y, with six turbulent 

forces F?,F?,Fy,FF,F?,Fy. 

We have solved (|3.23p ~ (|3.25|) numerically, using as inputs F u and F w obtained from 
our full simulations, i.e. the F modes shown in figure [TCJ The resulting solutions are 
shown in figure 1271 For comparison, we reproduce from figure [TU] the mean velocity 
fields, in Fourier representation, from our full simulations (DNS). The ODE solutions are 
virtually indistinguishable from the mean fields from DNS. Only in the sine component 
of U can the ODE solutions be distiguished (and only very slightly) from the DNS 
results. From the profiles in figure \%7\ the full mean fields could be constructed as in 
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Figure 27. Comparison between mean velocities (in Fourier representation) from full DNS and 
ODE models. Top row: Curves show U (blue, solid) and W (red, dashed) from DNS. Dots 
show solution to the full ODE model (|3.23|> . essentially indistinguishable from the solid curves. 
Bottom row: Curves again show U (blue, solid) and W (red, dashed) from DNS. Dots show 
solution to simplified ODE model (|3.26|l . The agreement with DNS is very good, though there 
are small differences particularly in the U component. 



figures [T2l and [T3"l Thus, while the ODE model requires input of the Reynolds-stress force 
terms, F u and F w , it demonstrates the simplicity of the force balance responsible for 
generating the patterned flow when viewed in the Fourier representation. Considering 
higher harmonics would be straightforward, but would serve little purpose. 

We can go in the other direction and attempt to simplify system (I3.23[) . The approxi- 
mate equalities w — Fq , F^ ss — F^ [see equation (|3. 12[1 ] . necessary for F to vanish 
at the center of the laminar region, can be imposed exactly, reducing the number of tur- 
bulent forcing input functions to four. The terms arising from the advective forces can be 
reduced by making approximations justified from figures [20l and [25l The nonlinear terms 
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m (|3.23a|) and (I3.23d|) can be neglected. The advective terms in (|3.236j) and (|3.23e[) can 
be approximated by ~k[3yU s and —k(3yW s . Making these approximations, we obtain: 



1 Ug + F% (3.26a) 



Re 



Q = -kPyU s + ^-U' ( !-FU (3.266) 
He 

= -V.{U' Q + a) + k(W + Py)U c + -^-U'J + Ff (3.26c) 

He 

= — WZ + F^ (3.26d) 
Re 

® = -kf3yW s + ^-W' c ! -F™ (3.26e) 
tie 

= -VM + 0) + k(W + (3y)W c + ±-W'J + (3.26/) 

He 

The solutions to this simplified ODE model are also presented in figure [27] There is 
quite good agreement with full DNS results, thus demonstrating that the dominant force 
balance is captured by this very simple system of ODEs. We stress that the only nonlin- 
earities in this model are in equations (|3.26c|) and fl3.26j| ). This reflects the complexity 
of the dynamics in the turbulent-laminar boundaries regions (and the simplicity of the 
dynamics in the centre of the turbulent and laminar regions.) 

From the simplified ODE model we can obtain the approximate equation satisfied at 
the centre of the laminar region by adding equations (|3.26a|) and (|3.26 6[) : 

k(3yU s = ^-(U Q + U c )". (3.27) 
Re 

This is a restatement in terms of Fourier components of the balance described by equation 
423) and figured 



4. Discussion 

We have presented an analysis of a particular turbulent-laminar pattern obtained in 
simulations of large-aspect-ratio plane Couette flow. We have focused on a single example 
so as to understand in quantitative detail the structure of these unusual flows. The 
key findings obtained in our study are as follows. First we find that in the (quasi-) 
laminar flow region the velocity profiles are not simply those of linear Couette flow. 
Instead a non-trivial flow is maintained in the laminar regions by a balance between 
viscous diffusion and nonlinear advection. Next we have considered the symmetries 
of the flow. When the pattern forms, the time-averaged flow breaks the translation 
symmetry but not centrosymmetry. The patterned state is centrosymmetric about the 
centre of the laminar region and about the center of the turbulent region. Next we 
have considered a spatial Fourier decomposition of the mean flow in the direction of the 
pattern wavevector. From this we find that the lateral structure of the pattern is almost 
completely harmonic, i.e. composed of a constant and single harmonic. Thus the pattern 
description can be reduced to just three cross-channel functions for each field, in that 
\J(x,y,z) w Uo(y) + U c (y) cos(fcz) + U s (y) sin(fcz). The absence of higher harmonics 
suggests that the pattern is near the threshold, in some sense, of a linear instability of 
a uniform turbulent state. Such an instability would be governed by a linear equation 
with coefficients which are constant in z, whose solutions are necessarily trigonometric 
in z. 
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From our analysis of the turbulent-laminar pattern, in particular its Fourier decompo- 
sition, we derive a model which reproduces the patterned flow. The model is derived from 
the averaged Navier-Stokes equations with the following assumptions. The crucial as- 
sumption, which is strongly supported by our numerical computations, is that the mean 
flow can be expressed in terms of just three horizontal modes. Effectively the generation 
of higher harmonics via nonlinear terms in the Navier-Stokes equations is negligible in 
the mean flow. The model is then further simplified because viscous diffusion is dom- 
inated by cross-channel diffusion - the standard boundary-layer approximation - and 
because pressure variation is negligible along the pattern wavevector. We take as input 
to the model the Reynolds-stress forces measured from computations. Assuming that 
the Reynolds stresses exactly vanish in the centre of the laminar regions, the number of 
inputs to the model is just four cross-channel functions. The result is a system of six 
simple ordinary differential equations which depend on four forcing functions. The model 
equations accurately reproduce the mean flow from full direct numerical simulations. 

A number of other researchers have attempted to reduce the description of turbulent 
or transitional plane Couette flow by various means. At these low Reynolds numbers, 
there is no doubt that fully resolved direct numerical simulation is feasible and gives 
accurate results. The purpose of formulating a reduced description is therefore to yield 
understanding. We now comment on the differences between the approaches used by 
other authors and our reduction. 

In parallel with their experiments, Prigent et al. (|2002l ; 120031 ) considered a pair of cou- 
pled Ginzburg-Landau (GL) equations with additive noise as a model for the transition 
from uniform turbulence to turbulent-laminar banded patterns via noisy (intermittent) 
patterns. These equations describe the variation in time and spanwise coordinate of the 
amplitudes A of two sets of laminar bands at opposite tilt angles. These laminar bands 
modulate the uniform turbulence in competition with one another. Each equation sepa- 
rately has one reflection symmetry which corresponds physically to the centrosymmetry 
K yz (equation (|3.6p ) of a banded pattern. The coupled GL equations possess a second 
reflection symmetry, corresponding physically to a spanwise reflection, which takes the 
amplitude A + to A~ and vice versa. By design, this symmetry is not present in our nu- 
merical computations. Prigent et al. used their experimental results to fit the parameters 
in the GL equation and then compared simulations of the equations with experimental 
results. Steady patterns in the resulting GL equations have only one non-zero amplitude 
and this amplitude possesses the reflection symmetry corresponding to n yz . Hence, the 
steady patterns in these simulations have exactly the symmetries of the patterns we have 
considered. 

An important class of models aims at reproducing dynamics of streamwise vortices and 
streaks in plane Couette turbulence by using a small number of ordinary differential equa- 
tions (ODEs). These equ ations de s cribe the time-evolution of amplitudes of modes with 
fixed spatial dependence. IWaleffd (jl997f) . guid ed by the discovery of the self-sustaining 
process (SSP) in direct numerical simulations ( Hamilton et q?J . Il995l ). derived a system 
of eight ODEs, whose variables represent amplitudes of the key ingredients of the SSP, 
namely longitudinal vor tices, streaks, and streak w aviness. This model was later also 
studie d and extended bv lDauchot fc Vioujardl ( 2000h and bv lMoehlis. Faisst fe Eckhardtl 
( 2004h . 

Two other Galerkin projection procedures have been used to derive ODE models. 
The most energetic streamwise-independent modes in a principle orthogonal composi- 
tion has been used as a basis for a 13-equation model ( Moehlis. Smith fe Holmesl l2002f) 
exhibiting heteroclinic cycles; when streamwise-depende nt modes are added, the result- 
ing 31-equation model ([Smith. Moehlis fe Holmesl . l2f)05h reproduces elements of the SSP 
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cycle . Eckhardt and co-workers ( Schrniegel fc EckhardtL 1997 : Eckhardt k. MersmannL 



1999) have proposed a Fourier space truncation of the Navier-Stokes equations in all three 
spatial directions leading to a 19-equation model. They calculated turbulent lifetimes 
and saddle-node bifurcations giving rise to new steady states in this model. 

Manneville and co-workers I Manneville fc LocherL [2000; Lag ha &: Mannevillel [2006) 



have proposed a drastic Galerkin truncation in the cross-channel direction y, retaining one 
or two trigonometric (for free-slip boundary conditions) or polynomial (for rigid boundary 
conditions) basis functions, but fully resolving both lateral directions. Simulating the 
resulting PDEs using a Fourier basis, they have been able to study phenomena such as 
the statistics of lifetimes of turbulent spots in domains with very large lateral dimensions. 

The reduction we have presented differs from the aforementioned studies in several 
respects. Most importantly, we do not describe any time-dependent behaviour. We 
consider here neither turbulent-laminar patterns which are themselves dynamic (as in 
Prigent et al.), nor do we consider the dynamics of streaks and vortices within the tur- 
bulence, nor do we consider the transient dynamics of turbulence. Instead we focus 
on the spatially periodic mean flow of steady turbulent-laminar patterns. While the 
turbulent portions of patterns are dynamic, containing streaks and streamwise vortices, 
these are on a fine scale relative to spatial scales of interest here. Our model description 
follows directly from an analysis of full numerical simulations (not from any a priori 
assumptions, physical or phenomenological), that show that all averaged velocity com- 
ponents and forces, including the Reynolds stress force, are almost exactly trigonometric 
in the direction of the pattern wavevector. It follows directly that the steady Reynolds- 
averaged Navier-Stokes equations can be reduced to 6 ODEs for cross-channel profiles of 
the Fourier modes. 

One of the more significant aspects of this work is the consideration of the force balance 
in just the laminar region. This balance is expressed by simple equations either in 
physical space, equation (|3.4|) . or in Fourier space, equation (|3.27|) . These equations are 
particularly interesting because they do not contain the Reynolds stresses, as these are 
negligable in the laminar region, and hence their implications can be understood without 
the need for closure assumptions. 

As noted in §3.11 equation (|3.4p implies that a non-zero tilt angle is necessary to 
maintain the S-shaped U profile in the laminar region. If the patterns were not tilted, the 
flow would necessarily be laminar Couette flow in the centre of the laminar regions where 
the turbulence vanishes. We can also derive implications for the relationship between 
Reynolds number, tilt angle and wavelength of the patterns from equation (|3.4j) . which 
we rewrite as: 

Re sin0 = (Uq + U c )" 

X 2n yU s 1 ■ ' 



Except where y U s ~ 0, the function on the right-hand-side is indeed approximately 
constant across the channel, between about 2.8 and 3.6. The value of Re sin#/A used in 
our simulations is 350 sin(24°)/40 = 3.56. 

We may obtain a qualitative understanding of this constant as follows; see figure [2H 
(left). Observe that in the center of the laminar region, the functional form of U = Uq+U c 
is like sin(7r?/). Hence its second y derivative can be approximated by multiplication by 
— 7r 2 , or equivalently (Uq + U c )"/(— it 2 ) w Uq + U c . We also find that the odd function 
—2 y U s is close to Uq + U c and is in fact almost indistinguishable from (Uq + U c )" /(— 7r 2 ). 
This implies that the right-hand-side of equation (|4.1| is nearly constant across the 
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Figure 28. Left: Comparison of profiles of Uq + U c (solid), (Uo + U c )"/(— n 2 ) (dashed) and 
—2yU s (points). The profile of (Uo + U c )"/(— tt 2 ) is close to Uo + U c , in accordance with 
the approximate functional form (Uq + U c ) ~ sin(7ry). Note that (Co + U c )" / (~tv 2 ) is almost 
indistinguishable from —2yU s , showing that (Uo + U c )"/(2n y U s ) is near ir over the entire 
channel. 

Right: Plot of Re sin#/A as a function of Re for the experimentally observed patterns of Prigent 
et al. ((2003, 200$. The open triangle shows Re sin(9/A = 3.56 for the case studied numerically 
in this paper. 



channel and equal to n, leading to: 

Re sin 6? . , „. 

— ^— a 7T (4.2) 

We believe that equation (|4.2p provides a good first approximation for the relation- 
ship between Re, A, and 9. Figure |2"51 (right) shows a plot of Re sin#/A as a function 
of Re from the experimental data of Prigent et al. (|2003r ). It can be seen that this 
combination of quantities is approximately constant with a value near 7r. The range of 
values of the individua l factor s Re, 9, and A can be seen in table [3] In prior studies 
(|Barklev fc Tuckerman! 12005 Jib , we have studied a large range of Reynolds numbers 
and tilt angles in a domain of length L z = 120. In this domain, the wavelength of a 
periodic pattern is less constrained, though it must be a divisor of 120. Figure |2T)1 shows 
the observed states as a function of Re and 9. Equation (|4.2p captures the correct order 
of magnitude of Re sin(9/A; specifically 1.8 < Re sinfl/A < 5. Moreover, in figure [29] one 
sees that for fixed Re, A increases with increasing 9, as (|4.2[) predicts. 

Equation (|4.2p does not hold in detail, however. Most notably, figure [25] shows that 
when Re is decreased at fixed 9, the wavelength A increases rather than decreases as one 
would expect from (|4.2p . We believe that the force balance (|3.4p holds for all patterns 
which possess a laminar region free of turbulence, but that the additional approximations 
made in deriving the simple relationship (|4.2p do not hold over the full range of conditions 
considered in figure [2"9"1 In particular, the right-hand-side of (|4.ip depends implicitly on 
Re, 9, and A via the dependence of Uq, U c , and U s on these quantities. The approximate 
functional relationships between Uq, U c and U s that we have observed in our simulations 
and on which we have relied in deriving (|4.2p may not hold for other parameter values. 
Finer adjustments must come from another mechanism. 

The main issue not addressed in our study is closure. We have not attempted to relate 
the forcing of the mean flow due to Reynolds stresses back to the mean flow itself. In 
the future we will report on studies employing closure models. 



We thank F. Daviaud, O. Dauchot, P. Le Gal, P. Manneville and A. Prigent for helpful 
comments. The simulations analyzed in this work were performed on the IBM Power 
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Figure 29. Patterns as a function of Reynolds number Re and 8 in a computational domain 
of siz6 Lx x Ly x L z — (4/sin(?) x 2 x 120. Turbulent-laminar patterns with wavelength A = 40 
(x), A = 60 (•), A = 120 (*). Uniform turbulence (■), intermittent turbulence (G3), laminar 
Couette flow (□). Wavelengths in computations are constrained to be divisors of 1 20. Numbers 
are wavelengths of experimentally observed patterns of Prigent et al. (2003, 2005). 
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Table 3. Turbulent-laminar banded patterns in plane Couette (PC), Taylor-Couette (TC) , 
rotor-stator (R S ), and plane Poiseuille (PP) flo w. Parameters reported in lPrigent et al\ ((2003) ; 
Cros fe Le Gall (|2002l );]Tsukaha ra et alj (|2005P ) are converted to a uniform Reynolds number 
based on the average shear and half-gap, as described in the Appendix. The two columns 
correspond to the values at the minimum and maximum Reynolds number reported. 
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Appendix. Turbulent-laminar bands in other shear flows 

Turbulent-laminar banded patterns have been observed in a number of shear flows: 
plane Couette (PC) flow, Taylor-Couette (TC) flow, rotor-stator (RS) flow (torsional 
Couette flow; the flow between differentially rotating disks) and plane Poiseuille (PP) 
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flow (channel flow). Comparisons between these flows are impeded by the fact that 
different conventions are used to non-dimensionalise each of them. 

In order to compare their observations in Taylor-Cou ette flow with those in plane 
Couette flow, iPrigent, Gregoire. Chate fc DauchotJ (|2003l) generalise the Reynolds num- 
ber used in plane Couette flow U = y/h by considering it as based on the shear and the 
half- gap: 



RePC = (Shear PC ) (half-gap) 2 



U 2 

h 



Uh 



(Al) 



V V V 

For flows whose shear is not constant, the average shear is used. We also convert stream- 
wise and spanwise wavelengths to total wavelength and angle of the pattern wavevector 
via 

As 



tan(0) = 



*span 



A 2 = A 



span 



cos(0) 



(A 2) 



Table [3] presents the Reynolds numbers, wavelengths, and angles for which turbulent- 
laminar patterns have been observed experimentally or numeric ally. The sub s ection s 
which follow explain h ow Table El was obtaine d from the data in IPrigent et all ([20031 ) ; 
Cros fc Le Gal |2002h : iTsukahara et all (|2005l ). 



A.l. Taylor- Couette flow 

For Taylor-Couette flow between differentially rotating cylinders, the azimuthal and axial 
directions correspond to the streamwise and spanwise directions of plane Couette flow. 
For cylinders of radius and r a , rotating at angular velocities u>i and lo with 2h = r —Ti 
and r] = Ti/r , the shear averaged over the gap is 



nuii - rjr u; 



(1 + V)h 



(Shcar TC ) = 
leading to the Reynolds number: 

TiUi - i]r uj h 2 Rei - Re Q 



Re 



TC 



+ v 



(A3) 



(A 4) 



where the last approximate equality corresponds to exact counter-rotation (u) a = —u>i) 
and the narrow gap limit (77 — » 1), and Ri, R a are the conventionally defined inner 
and outer Reynolds numbers, e.g. Rei = 2hriLUi/is. The wavelengths and Reynolds 
n umbers observed in T aylor-Couette and plane Couette flow are compared in Figure 5 
of IPrigent efal] (|2003l ). 



A. 2. Torsional Couette flow 

The laminar profile for torsional Couette flow between a rotating and a stationary disk 
(rotor-stator flow) is 

ujr z 



and the Reynolds number based on axial shear and half-gap is 

,2 

Av 

For m spirals, the azimuthal wavelength in units of the half-gap is 



uir h 2 Lorh 
h Av 



A 



RS 



27IT 

mh/2 



Airr 
mh 



(A 5) 
(A 6) 
(A 7) 
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Turbulent spiral patterns which are rath er regular occur for a range of angular velocities 
and radii. In their figures 12, 16 and 18. ICros fc Le Call (|2002h focus particularly on the 
radius and gap: 

r = 0.8 x 140 mm = 11.2 cm h = 0.22 cm (A 8) 

The highest and lowest rotation rates for which turbulent spirals are seen are 

(J = 68rev/min = 7.12rad/sec with m = 6 (A 9a) 

ui = 47rev/min = 4.92rad/sec with m = 9 (A96) 

Substituting fO|) - l|X9|) and the viscosity v = 10~ 2 cm 2 /sec f water into (|A~6l) - ([A~7f 
leads to the values shown in Tabled The pitch angle of the spirals remains approximately 
constant at 19°. We use (|A 2\ to calculate A span and A 2 , neglecting the variation in radius. 



A. 3. Plane Poiseuille flow 

Figure 14 of iTsukahara. Seki. Kawamura fc Tochiol (j2005h shows a visualisation from a 
direct numerical simulation of plane Poiseuille (PP) flow in a channel with domain and 
Reynolds number 

u 5 

^stream x L y x L span = 51.2 8 x 2 8 x 22.5 8 Re c = — = 1430 (A 10) 

where u c is the centerline velocity. The domain contains a single wavelength of an oblique 
turbulent-laminar banded pattern oriented at 8 = 24° to the streamwise direction. (Both 
the wavelengt h and the angle are dictated by the computational domain.) Following 



Waleffd ([20031 ). we view the Poiseuille profile in the half-channel [—8, 0], over which the 



shear has one sign, as comparable to the Couette profile in the channel [—h, h], and thus 
take 8/2 as the unit of length, rather than 8. The shear is obtained by averaging over 
H,0]: 



pp. / du \ 



(Shear"*) = = -f (All) 

For the Reynolds number based on the average shear and half-gap, we obtain 
RcPP _ (Shear pp ) (half-gap) 2 _ l -f S -^ _ u c 8 _ Re c _ 1430 _ ^ g 
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